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Abstract 

We present the Fourier Transform of a continuous gravitational wave. We have 
analysed the data set for one day observation time and our analysis is applicable for 
arbitrary location of detector and source. We have taken into account the effects arising 
due to rotational as well as orbital motions of the earth. 



1 Introduction 

The first generation of long-baseline laser interferometers and ultra cryogenic bar detec- 
tors will start collecting data very soon. The network of detectors is expected not only to 
confirm the existence of gravitational waves (GW) but will also provide the information 
about the structure and the dynamics of its source. At the present stage, the data analysis 
depends largely upon the study of the expected characteristic of its potential sources and 
the waveforms. Majority of the experimental searches is focused on the detection of burst 
and chirp signals. However, the interest in the data analysis for continuous gravitational 
wave (CGW) signals is growing. A prime example of sources of this type is a spinning 
neutron star. Many research groups around the globe are working extensively on the data 
analysis for spinning neutron stars (Jaranwoski, et.al 1998, 1999, 2000; Brady et. al. 1998, 
2000; Krolak 1999). 
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The detection of GW signals in the noisy output of the detectors has its own problems, 
not the least of which is the sheer volume of data analysis. Bar detectors have essentially 
the same problems as interferometers in reference to CGW sources. Each detector pro- 
duces a single data stream that may contain many kinds of signals. Detectors don't point, 
but rather sweep their broad quadrupolar beam pattern across the sky as the earth moves. 
So possible sources could be anywhere on the sky and accordingly the data analysis algo- 
rithms need to accommodate signals from any arbitrary location of its source. 

In this and the subsequent paper we present analysis of Fourier Transform (FT) of the 
output data of a ground based laser interferometer. The output data has broad band noise 
and the signal is to be extracted out of it. For this, one has to enhance signal-to-noise 
ratio (SNR). This is achieved by analyzing long observation time data as SNR is directly 
proportional to the square root of observation time (V^)- However, in a data for long 
duration, the monochromatic signal gets Doppler modulated due to (i) the orbital motions 
of the earth around the sun and (ii) the spin of the earth. The frequency modulation (FM) 
will spread the signal in a very large number of bins depending on the source location and 
the frequency. In addition to this, there is amplitude modulation (AM). As we will see 
in the sequel the amplitude of the detector output consists of simple harmonic terms with 
frequencies Wrot and 2wj.ot where, Wrot stands for angular rotational frequency of the earth. 
Accordingly, the AM results in splitting of the FT into frequencies ±Wrot and ±2wrot- 

In the next section we present the noise free response of the laser interferometric de- 
tector and obtain the explicit beam pattern functions. In section 3 we discuss the Doppler 
effect and obtain Fourier transform (FT) of the FM signal for arbitrary source and detector 
locations taking into account the earth's rotational motion about its axis and its revolution 
around the sun. In section 4 the FT of the Doppler modulated complete response of the 
detector is obtained. Section 5 is the conclusion of the paper. 

2 The noise free response of detector: Beam pattern and 
Amplitude modulation 

Let a plane GW fall on a laser interferometer and produce changes in the arms of the 
detector. In order to express these changes quantitatively we have to specify the wave and 
the detector. Let XYZ and xyz represent respective frames characterising the wave and 
the detector. We assume the direction of propagation of the wave to be the Z axis and the 
vertical at the place of the detector to be the z axis. The difference of the changes SI in the 
arm lengths of the detector may be given via 

R(t) = - = - sin 20 [(A^A^ - A^A^) h+ + (A^A^ + A^A^) h^] (1) 
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where is the normal length of the arms of the detector and 2f2 expresses the angle be- 
tween them (Schutz and Tinto, 1987). The matrix (A^^-^ represents the transformation 
expressing the rotations to bring the wave frame (X, Y, Z) to the detector frame (x, y, z). 
The direction of the source may be expressed in any of the coordinates employed in Spher- 
ical Astronomy. However, we find it convenient to define it in Solar System Barycentre 
(SSB) frame, (X', Y\ Z'). This SSB frame is nothing but astronomer's ecliptic coordinate 
system. Let and denote the celestial colatitude and celestial longitude of the source. 
These coordinates are related to the right ascension, a and the declination, 5 of the source 
via 



cos 9 = sin S cos e — cos S sin e sin a 
sin 9 cos (p = cos 6 cos d 
sin 9 sin 6 = sin S sin e + cos S cos e sin a 



(2) 



where e represents obliquity of the ecliptic. We choose x axis as the bisector of the angle 
between the arms of the detector. At this stage the orientation of the detector in the hori- 
zontal plane is arbitrary. It is assigned with the help of the angle 7 which the x axis makes 
with the local meridian. The location of the detector on the earth is characterised by the an- 
gles, a - colatitude and P - the local sidereal time, expressed in radians. The transformation 
matrix (A]^^ may be expressed as 



A = DCB 
where 

B : rotation required to bring XYZ to X'Y'Z' 
C : rotation required to bring X'Y'Z' to x'y'z' 
D : rotation required to bring x'y'z' to xyz 



(3) 



Here x'y'z' represents the frame associated with the earth. The Euler angles defining the 
corresponding rotation matrices are given via 



B 
C 
D 



(0,e,0) 
(a,/5 + 7r/2,7-7r/2) J 



(4) 



where is a measure of the polarisation of the wave (Goldstein, 1980). 
Let us write Eq. ([1|) as 



m = r 



sm2n [F+h+ + Fx/i> 



(5) 



The functions Fj^ and Fx involve the angles 9, (f), ip, e, a, P, 7 and express the effect of 
the interaction of the wave and the detector. These are called antenna or beam patterns. 
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The explicit functional dependences of F+ and Fx are given by Jotania and Dhurandhar 
(1994). It is pointed out that there are few minor errors in these expressions and these are 
later corrected by the authors in Jotania (1994). These functions appear complicated but 
may be written in simpler form by introducing following abbreviations. 

U — cos a cos (3 cos 7 — sin /3 sin 7 , 

V — — cos a cos /3 sin 7 — sin /3 cos 7 , 
X — cos a sin /3 cos 7 + cos /3 sin 7 , | 

Y — — cos a sin /3 sin 7 + cos P cos 7 



L = cos ip cos — cos 9 sin sin -0 , 

M = cos 'ip sin (i) + cos9 cos ci sin ■(/' , 

N = — sin '0 cos — cos^^ sin cos ■0 , 

P = — sin -0 sin + cos^ cos cos '0 , 

Q — sin ^ sin , i? = sin 9 cos , 



(7) 



A 
B 
C 
D 
E 



2XY cos^ e — sin^ e sin^ a sin 27 + sin 2e(X sin a sin 7 — F sin a cos 7) , 
2XY sin^ e — cos^ e sin^ a sin 27 — sin 2e{X sin a sin 7 — y sin a cos 7) , 
cos e{YU + XV) + sin e{U sin a sin 7 — sin a cos 7) , 
— sin €{YU + XV) + cos e(C/ sin a sin 7 — y sin a cos 7] , 

— 2Xy cos e sin e — cos e sin e sin^ a sin 27 + cos 2e{X sin a sin 7 — y sin a cos 7) 



>(8) 



After straight-forward substitutions one obtains 



F+{t) = - [2(L2 - M^)UV + {N^ - P^)A + - R'^)b] + {LN - MP)C 

+ {LQ + MR)D + {NQ + PR)E , (9) 

Fx (t) = 2LMUV + NPA - sin^ 9 sin 20 + {LP + MN)C 

+ (Mg - LR)D + (Pg - iVi?)F (10) 



The compactification obtained here arises due to the fact that above introduced abbrevia- 
tions find places in the transformation matrices as follows 
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B 



L 

M 
sin 6 sin ip 



N 
P 

sin 9 cos ip 



Q 

-R 
cos 9 



(11) 



1 
COS e sin e 
— sin e cos e 



(12) 



U F sin a cos P 

X Y sin a sin (3 j , 

sin a cos 7 sin a sin 7 cos a 



(13) 



After algebraic manipulation Eqs. and ( [T^ ) may be expressed as 



F+ (t) = Fi^ cos 2/3 + sin 2/3 + Fg^ cos /3 + F^^ sin /3 + Fg^ ; 



Fx (t) = Fi, cos 2/? + F2, sin 2/? + F3, cos p + F4, sin /? + F5 



(14) 
(15) 



where Fj^ and Fj^ (2 = 1, 2, 3, 4, 5) are time independent expressions given via 

Fi^ = -2Gcosacos27 + ^^^^f^(cos2a + 1) , ^ 

F2^ = if COS a COS 27 + G sin 27(cos^ a + 1) , 

Fs^ = / sin a cos 27 + J sin 2a sin 27 , 

-P4+ = 2J sin a cos 27 — I sin 2a sin 27 , 



(16) 



G = 1[{LQ + M R) sin e-{LN -MP) cose], 

H = \ [{N^ - P2) cos2 e - (L2 - M^) + (Q2 _ ^^2) ^^^2 ^ _ ^j^q ^ gij^ 2e] , 

/ = i[(Q2 _ gij^2e - (iV2 - p2) sin2e -2(iVg + Pi?) cos 2e], 

J = i [(LiV- MP) sin e + (LQ + MP) cose] 

Let us note that Fj^ is related to Fj_^ via 
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I = 



1,2,3,4,5 



(18) 



This symmetry is representative of the quadrupolar nature of the detector and the wave. 
The detectors at different orientations will record different amplitudes in their responses. 
The explicit beam pattern functions may be computed easily for any instant of time. Due to 
the symmetries involved in and Fx it is sufficient to evaluate either of the beam patterns. 

The amplitude modulation of the received signal is a direct consequence of the non- 
uniformity of the sensitivity pattern. As remarked earlier they are fairly complicated func- 
tion of their arguments. Equations ( p^ and ( p3| ) reveal that the monochromatic signal 
frequency will split, due to AM, into five frequencies. This results in the distribution of 
energy in various frequencies and consequent reduction of the amplitude of the signal. The 
periodicity of the beam pattern and Fx with a period equal to one sidereal day is due to 
the diurnal motion of the Earth. 

3 Doppler shift and Frequency modulation 

The frequency of a monochromatic signal will be Doppler shifted due to the translatory 
motion of the detector acquired from the motions of the earth. Let us consider a CGW 
signal of constant frequency /„. The frequency /' received at the instant t by the detector 
is given by 



where n is the unit vector from the antenna to the source, v is the relative velocity of the 
source and the antenna, and c is the velocity of light. The unit vector n from the antenna to 
the source, in view of the fact that the distance of the source is very large compared to the 
average distance of the centre of the SSB frame and the detector, may be taken parallel to 
the unit vector drawn from the centre of the SSB frame to the source. Hence, 

n = (sin 6 cos (p , sin 6 sin , cos 6) (20) 

As V keeps on changing continuously both in its amplitude and direction /' is a continuous 
function of t. Further, since f ^ c we take 7o = 1. 

The radius vector r{t) in the SSB frame is given by (Kanti, et al. 1996) 

r(t) = [Rse cos{worbt) + Re sin a cos P , Rse sm{worbt) + Re sin a sin P cos e — 




(19) 



Re COS a sin e , Re sin a sin /? sin e + Re cos a cos e] ; 



(21) 
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(3 ^ (3„ + Wrott 



(22) 



where Re, Rse and Worb represent respectively the earth's radius, average distance between 
earth's centre from the origin of SSB frame and the orbital angular velocity of the earth. 
Here t represents the time in seconds elapsed from the instant the sun is at the Vernal 
Equinox and /3o is the local sidereal time at that instant. The Doppler shift is now given via 



f-fo 
fo 



v.n 



it) 



rt.n RseWorb 



sin 9 sin(0 — Worht) + 



sin a [sin ^{cos (3 cos e sin — cos sin /3} + cos (3 sin e cos 9] (23) 



The phase of the received signal is given by 



m) = 2ti t f{t')dt' 

Jo 



JO 



v.n 



it') 



dt' 



(24) 



Here we assume the initial phase of the wave to be zero. After straight-forward calculation 
we obtain 



R R 
t H — sin 9 cos 0'H sin alsin ^(sin /3 cos e sin + cos cos /3) + 

c c 



sin (3 sin e cos 9} 



Rx 



sin 9 cos ( 



Re 



sin a{sm 6* (sin /3o cos e sin + 



cos (f) cos f3o) + sin f3o sin e cos ^^}] 

27r/ot + Z COs{Worbt — (f)) +V Sm{Wrott) + Q COs{Wrott) — Ti- — Q 

27r/ot + Z cos{a^rot - 4>) + cos{^rot -S) -U- Q 



(25) 



where 



V 
Q 

u 

z 
n 



27: sin a (cos /3o(sin 9 cos e sin + cos 9 sin e) — sin Po sin ^ cos 0) , 
27r/o— sin Q;(sin /5o(sin ^ cos e sin + cos 9 sin e) + cos /3o sin ^ cos 0) , 



27Tfo^sm9, 
Z cos , 



(26) 
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tan" 



-1 £ 

Q 



<P' = Worbt - <P , 

Lrb = Worbt = a^roU a = Worb/Wrot ~ 1/365.26, 
^rot ^rott 

The two polarisation states of the signal can be taken as 

h+{t) cos[$(t)] 
h^{t)=K, sin[$(t)] 

ho+, ho^ are the time independent amplitude of h+{t), and /ix (t) respectively. 
To understand the nature of the FM let us consider the function 

h{t) = cos[<l>(t)] 

and analyse it for one day observation data. The FT is given via 



Hf) 



T 



<i Jo 

This splits into two terms as 



cos[$(t)]e ; T = one sidereal day = 86164 sec. 



Hf) 



(27) 



(28) 
(29) 



(30) 



(31) 



(32) 









' ! 




2Wrot JO 






I^. = 


' f 

2Wrot JO 




foTf_ 




frot 



di. 



2-K 



-i[iu++Z cos(ai,-<i>)+M cos(i-5)-n-Q] 



) ^ ^rot l^rott 



(33) 
(34) 
(35) 



Numerical result shows that I^^ oscillates very fast and contributes very little to h{f) 
Hence, hereafter, we drop from Eq. ( p2| ) and write v in place of z/_. Using the identity 



/=oo 



g±*Kcostf = Jo{±k) +2J2 i^M^f^) COS W 

1=1 



(36) 
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we obtain 



A- 



. /•27r 

Jo 



k=oo 



k=l 



2Wrot 



X 



Jo{Z) + 2 ^ Jk{Z)i^ cos k{ai - 0) 



m=l 



(37) 



where J stands for the Bessel function of the first kind. After performing the integration 
we get 



Kf) 



k=oo m=oo 



, 2^ E E e^Me-iv\-. 

^Wrot k=-oo m=-oo 



(38) 



A = 

B = 

C = 

V = 



n-Q 



ak+m 



{cos 2i/7rsin(2aA;7r — k(t) — m5) + 



Jk(Z)Jm(M) 

sin 2i/7r cos(2a/c7r — kef) — mS) 
sm{k4> + m5)} 

cos 2i/7r cos(2aA;7r — kcj) — mS) + sin 2i/7r sin(2aA;7r — kef) — m5) - 

cos{k(f) + m5) 



(39) 



The FT of the two polarisation states of the wave can now be written as 



ho^ [Kf)] 



k=oo m=oo 



"0+ 



2w. 



(40) 



rot m=—<x> 



h.if) 



iho^ [Kf) 



I'hn 



k=oo m=oo 



2m 



rot f^- 



(41) 



-oo m=— oo 



The FT of the FM signal contains the double series Bessel functions. The Bessel functions 
have contributions due to the rotational as well as the orbital motion of the earth. It is 
remarked that Jotania et al. (1996) have analysed FT of FM signal for one day observation 
time. They have taken specific detector as well as source location. They have also neglected 
the orbital motion. Our analysis generalizes their results. We may now compute h{f) 
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and may plot its behaviour. To achieve this we have made use of Mathematica software. We 
know that the value of Bessel function decreases rapidly as its order exceed the argument. 
Accordingly, one computes in practice only a finite number of terms in the above infinite 
series. Figure ([1]) represents such a plot for 



fo = 80 Hz, ho^ = ho, = l 

a = n/3, /5o = 7r/4, 7 = 27r/5, 

6 = 7t/36, (p = it , il) = 7i/6. 



(42) 



with a resolution equal to 1/To = 1.16 x 10^^ Hz. For the present case we found it 
sufficient to evaluate the infinite series for k = —21900 to +21900 and m = —10 to +10. 
Figures @ and (^ represent the plot of the FT at resolution 10^^ Hz and 10^^ Hz. 
A careful look at these plots reveals that the resolution of Fig. ([T|) does not represent the 
details of the dominant peaks around fo, whereas, Fig. (Q) does not give any new behaviour 
as compared to Fig. @. Hence, we may say that a resolution of about 10^^ Hz is required to 
understand the correct behaviour of the FT for one day observation data. In this reference 
let us recall that the data analysis for Fast Fourier Transform (FFT) limits the resolution 
to I/Tq. However, the detector output may provide us higher resolution. Thus the semi- 
analytical analysis presented here may provide more information as compared to FFT. 



4 Fourier transform of the complete response 

The complete response R{t), in view of Eqs. (^, (jT^, (jT3|), (^ and (|2^) may be written 
as 



R{t) = R+{t) + R^{t) ; 



(43) 



Fi^ cos 2(3 + sin 2(3 + F^^ cos (3 + sin (3 + Fg^ 
Fi^ cos 2/? + F2^ sin 2/5 + Fg^ cos/? + sin/5 + Fg^ 

Here, for simplicity, we have taken the angle between the arms of the detector equal to 7r/2 
i.e. n = Tx /A. Now the FT of the complete response may be expressed as 



R+{t) = ho 
Rx{t) = K 

re, for si 

Q = IT/ 



cos[<l>(t)], (44) 
sin[$(t)] (45) 



Rif) = RM) + RAf) 

Substituting /3 as given by ( p2| ) one obtains 



(46) 



e-'''^''iF,^+iF2^)e 



e-*/^o(F3^ + iF4^)e-^"''-°'* + e*^°(F3^ - iFi^)e''"^°'' + 2F^^ cos[$(t)] (47) 



-i2Wrott 



iWrott 
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and similar expression for i?x (t). Now it is straight-forward to obtain the expressions for 

i?+(/)andi?x(/) as 



''0+ 

2 



I rot; 



4+. 



h{f - frot) 



+ 2F, 



5+ 



(48) 



^x(/) 



''Ox 

2 



i2F, 



(49) 



Collecting our results the FT of the complete response of the detector for one day time 
integration will be 



1 r _ 

- e 

2 L 



- 2/,„t)]^ - ^^2+) - + 

/^(/ - /rot)] ^ [ V (^3+ - ^i^4+ ) - /^o. (i^4. + ^i^3 J 



+ 



+ 



+ 



(50) 



This shows that due to AM every Doppler modulated FM signal will split in four ad- 
ditional lines at / ± 2frot and / ± f^ou where frot is the rotational frequency of Earth 
Urot ~ 1.16 X 10-5 Hz). 



We have plotted in Fig. (Q) the power spectrum of the noise free complete response 
of the signal for its various parameters as given by (p2|). The contribution in the power 
spectrum of the modulation at frequencies / + 2/rot, / — 2jrou f + frot and / — frot and 
/ are represented respectively in Figs. (Q), @, (|8|) and @. It is observed that the most 
of the power will be at / + 2/rof and least power will be in / — frot- 
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5 Conclusion 

In this paper we have considered the effect of earth's motion on the response of the detector 
through FT analysis. It can be easily inferred from Eqs. ( pO| ) and ( pO| , |4T| ) that the splitting 
of frequencies (i) in AM arises explicitly due to rotational motion and (ii) in FM arises due 
to rotational as well as orbital motion of the earth. In view of the fact that the data output at 
the detector is available in discrete form, the analytical FT is not very convenient and one 
normally employs the popular FFT. However, FFT has resolution limited to 1 /To. Further, 
it is important to understand for how much time one can ignore the frequency shift arising 
due to Doppler effect. In fact, Schutz (1991) has demonstrated that these effects due to 
rotational motion are important after the time given by 



This means that for GW signal for frequency 80 Hz one has to take into account these 
effects after data time ^ 4 hours. The analytical FT studied in this paper leads to following 
inferences: 

(i) FFT for one day observation data will not provide sufficient resolution as to represent 
the correct picture of the frequency splitting. 

(ii) The adequate resolution required for one day observation ~ 10~^ Hz. 

(iii) The frequency split due to FM for frequency fo = 80 Hz and source at {9, 0) = 
(tt /36, vr) is ^ 2 X 10~^ Hz and due to AM is ^ 4.64 x 10"^ Hz. 

(iv) The drop in amplitude due to FM alone is about 56% 

(v) The drop in amplitude due to AM alone is about 18%. 

(vi) The drop in amplitude for the complete response is about 74%. 

(vii) The maximum power due to AM is associated with fo + 2frot- 

It is remarked that the drop of the amplitude in complete response is severe both due to AM 
and FM as the relevant frequency range lies in the same region. 

We have presented the FT analysis assuming the phase of the GW to be zero at that 
instant t = 0. However, one may relax this condition and may obtain the results easily by 
taking into consideration the effects of change of the time origin. 




(51) 
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Figure 1: FT of a FM signal frequency, = 80 Hz, from a source located at (7r/36, tt) 
with a resolution of 1.16 x 10~^ Hz 
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Figure 2: FT of a FM signal frequency, = 80 Hz, from a source located at {n/36, n) 
with a resolution of 10"^ Hz 
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Figure 3: FT of a FM signal frequency, = 80 Hz, from a source located at {n/36, n) 
with a resolution of 10~^ Hz 
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Figure 4: Power spectrum of the complete response of a Doppler modulated signal fre- 
quency, fo — 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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Figure 5: Power spectrum of a Doppler modulated signal at frequencies / + 2/^0* of signal 
frequency, = 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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Figure 6: Power spectrum of a Doppler modulated signal at frequencies / — 2fj.ot of signal 
frequency, = 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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Figure 7: Power spectrum of a Doppler modulated signal at frequencies / + frot of signal 
frequency, = 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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Figure 8: Power spectrum of a Doppler modulated signal at frequencies / — frot of signal 
frequency, = 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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Figure 9: Power spectrum of a Doppler modulated signal at frequencies / of signal fre- 
quency, fo — 80 Hz, from a source located at (7r/36, tt) with a resolution of 10~^ Hz 
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